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The N = 4 super Yang-Mills plasma is studied in the regime of weak coupling. Collective 
excitations and collisional processes are discussed. Since the Keldysh-Schwinger approach is used, 
the collective excitations in both equilibrium and non-equilibrium plasma are under consideration. 
The dispersion equations of gluon, fermion, and scalar fields are written down and the self-energies, 
which enter the equations, are computed in the Hard Loop Approximation. The self-energies are 
discussed in the context of effective action which is also given. The gluon modes and fermion ones 
appear to be the same as those in the QCD plasma of gluons and massless quarks. The scalar 
modes are as of free relativistic massive particle. The binary collisional processes, which occur at 
the lowest nontrivial order of the coupling constant, are reviewed and then the transport properties 
^^ of the plasma are discussed. The J\f — 4 super Yang-Mills plasma is finally concluded to be very 

similar to the QCD plasma of gluons and light quarks. The differences mostly reflect different 
^H numbers of degrees of freedom in the two systems. 

C\l I. INTRODUCTION 

Supersymmetric models are considered as possible extensions of the Standard Model, see e. g. PQ, and supersymme- 

try is then assumed to be a symmetry of Nature at a sufficiently high energy scale. Experiments at the Large Hadron 

O ! Collider may soon show whether this is the case. This paper is devoted to a plasma system with the dynamics governed 

by the Af = 4 supersymmetric Yang-Mills theory [21 H]. In the models with an extended (M > 1) supersymmetry, 

the left- and right-handed fermions interact in the same way, in conflict with the Standard Model where the left- and 

right-handed matter particles are coupled differently. Consequently, the M = 4 super Yang-Mills is not treated as a 

£f~) serious candidate for a theory to describe the world of elementary particles. Nevertheless, the theory attracts a lot of 

attention because of its unique features. The J\f = 4 super Yang-Mills appears to be finite and thus it is conformally 

>0 invariant not only at the classical but at the quantum level as well. 

A great interest in the AA = 4 super Yang-Mills theory was stimulated by a discovery of the AdS/CFT duality of 
the five-dimensional gravity in the anti de Sitter geometry and the conformal field theories :.4J , for a review see [5] and 
the lecture notes [6] as an introduction. The duality offered a unique tool to study strongly coupled field theories. 
Since the gravitational constant and the coupling constant of dual conformal field theory are inversely proportional 
to each other, somc problems of strongly coupled field theories can be solved via weakly coupled gravity. In this way 
some intriguing features of strongly coupled systems driven by the M = 4 super Yang-Mills dynamics were revealed, 
see the reviews [7J [S], but relevance of the results for non-supersymmetric systems, which are of our actual interest, 
remains an open issue. In particular, one asks how properties of the N = 4 super Yang-Mills plasma (SYMP) are 
related to those of quark-gluon plasma (QGP) studied experimentally in relativistic heavy-ion collisions. While such 
a comparison is, in general, a difficult problem, some comparative analyses have been done in the domain of weak 
& coupling where perturbative methods are applicable [5rTH]. 

We undertook a task of systematic comparison of supersymmetric perturbative plasmas to their non-supersymmctric 
counterparts. We started with the J\f = 1 SUSY QED, analyzing first collective excitations of ultrarelativistic plasma 
which, in general, is out of equilibrium [15j . We computed the one-loop retarded sclf-cncrgics of photons, photinos, 
electrons and selectrons in the Hard Loop Approximation using the Keldysh-Schwinger formalism. The self-energies, 
which we also analyzed in the context of effective action, enter the dispersion equations of photons, photinos, electrons 
and selectrons, respectively. The collective modes of Af = 1 SUSY QED plasma appear to be essentially the same 
as those in ultrarelativistic electromagnetic plasma of photons, electrons and positrons. In particular, a spectrum of 
photino modes coincides with that of quasi-electrons. Therefore, independently whether photon modes are stable or 
unstable, there are no unstable photino excitations. The supersymmetry, which is obviously broken in the plasma 
medium, does not induce any instability in the photino seetor. 

In the subsequent paper [16] we discussed collisional charaeteristies of M = 1 SUSY QED plasma. For this purpose 
we computed cross seetions of all elementary processes which occur at the lowest non-trivial order of a = e 2 /4tt. We 
found that some processes, e. g. the Compton scattering on selectrons, are independent of momentum transfer. The 
processes are qualitatively different from those of usual electromagnetic interaetions dominated by small momentum 
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transfers. Further on we discussed collisional characteristics of equilibrium Af = 1 SUSY QED plasma, observing 
that parameters of ultrarelativistic plasmas are strongly constrained by dimensional arguments, as the temperature 
is the only dimensional quantity of equilibrium system. Then, transport coefficients like viscosity are proportional to 
appropriate powers of temperature and the coefficients characterizing diffcrent plasmas can differ only by numerical 
factors. So, we derived the energy loss and momentum broadening of a particle traversing the equilibrium plasma, 
which dcpcnd not only on the plasma temperature but on the energy of test particle as well. We found that the 
two quantities have very similar structure (in the limit of high energy of the test particle) even for very different 
elementary cross sections. Our findings presented in [TS1[TB] show that the plasmas of N = 1 SUSY QED and of QED 
are surprisingly similar to each other. In this paper we discuss properties of the Af = 4 super Yang-Mills plasma, 
analyzing both collective excitations and collisional characteristics of the system. 

Our main aim is to confront the weakly coupled plasma driven by Af = 4 super Yang-Mills with the perturbative 
quark-gluon plasma governed by QCD. We do not attempt to compare our results to those obtained in strong coupling 
regime using either the AdS/CFT duality or lattice QCD. Some plasma characteristics we discuss, e. g. the energy loss, 
are computed in both strongly and weakly coupled systems but it is rathcr unclear how to study collective excitations 
representing colored quasiparticles in the setting of AdS/CFT duality or lattice QCD. The paper [T7] demonstrates 
that even the definition of Debye screening mass, which has a very simple meaning in perturbative plasmas, is far 
not straightforward in strongly interacting systems. For these reasons we escape from discussing our results in the 
context of strong coupling. 

Our paper is organized as follows. In the next section, we discuss the Lagrangian of Af = 4 super Yang-Mills and 
the ficld content of the system under consideration. The vertexes of Af = 4 super Yang-Mills are collected in Appendbc 



|A| In Sec. III basic characteristics of SYMP such as energy density and Debye mass are discussed and compared to 
those of QGP. Then, we move to plasma collective excitations. The general dispersion equations of gauge bosons, 
fermions and scalars are written down in Sec. |IV| and the self-energies, which enter the equations, are obtained in 
the subsequent section. We apply here the Keldysh-Schwinger approach which allows one to study equilibrium and 
non-equilibrium systems. The free Green's functions of Keldysh-Schwinger formalism are given in Appendix|BJ Since 
we are interested in collective modes, the self-energies are obtained in the long wavelength limit corresponding to 
the Hard Loop Approximation. The effective action of the Hard Loop Approach is derived in Sec. |VI| and possible 



structures of self-energies are considered in this context. In Sec. VII we present a qualitative discussion of collective 
modes in SYMP. Sec. |VIlJ] is devoted to collisional characteristics of the plasma - elementary processes and transport 
coefficients are briefly discussed here. Finally, we conclude our study in Sec. |IX| 



As we have intended to make our paper complete and self-contained, there is inevitably some repetition of the 



content of our previous publications |15l I16j. mostly in Sees. VII VIII Throughout the paper we use the natural 



system of units with c = % = ks = 1; our choice of the signature of the metric tensor is (H ). 

II. Af = 4 SUPER YANG-MILLS THEORY 

We start our considerations with a discussion of the Lagrangian of Af = 4 super Yang-Mills theory [3J [5] . We follow 
here the presentation given in [15] . 

The gauge group is assumed to be SU(Y C ) and every field of the Af — 4 super Yang-Mills theory belongs to its 
adjoint representation. The field content of the theory, which is summarized in Table ||J is the following. There are 
gauge bosons (gluons) deseribed by the veetor field A® with a, b, c, • • • = 1,2,... N% — 1- There are four Majorana 
fermions represented by the Weyl spinors A Q with a = 1, 2 which can be combined in the Dirac bispinors as 

*=(T), * = (A Q ,A«), (1) 

where A^ = [A Q ]^ with f denoting Hermitian conjugation. To numerate the Majorana fermions we use the indices 
i, j = 1,2,3,4 and the corresponding bispinor is denoted as ^i. Finally, there are six real scalar fields which are 
assembled in the multiplet $ = (Xi, Yy, X%, Y2, Y3, Y3). The components of $ are cither denoted as X p for scalars, 
and Y p for pseudoscalars, with p,q = 1,2,3 or as <&a with A, B = 1, 2, ... 6. 
The Lagrangian density of Af = 4 super Yang-Mills theory can be written as 

C = ~\f» v F; v + !*?(W + \{D^ A ) a {D^ A ) a (2) 



f a bef cde$a A $*,$*$* - z|f bc (*?«£. X b p ^ + i^fl^Y**'. 



where Fg v = d»A v a - d" 'A» + gf abc A£A" c and the covariant derivatives equal {lp%) a = {$5 ab + gf abc 4^ b and 
(D»$) a = D% b $ b = {d»S ab + gf abc A»)<S> b ; g is the coupling constant; f abc are the structure constants of SU(iV c ) 
group; the 4x4 matrices a p , f3 p satisfy the relations 

{aW} = -2<^, {P,IF} = -2S", [a", £«] = (), (3) 

and their explicit form can be chosen as 

^(17)- " 2 -(-°.°> "'-(Tta). < 5 > 

where the 2x2 Pauli matrices read 

u A 0\ , _ /0 1\ 2 _ /0 -A 3 /'i f) 



As seen, the matrices a p ,f3 p are antiHermitian: (a p )^ = — a p , (/? p )^ = — /3 P . The vertexes of j\f = 4 super Yang-Mills, 
which can be inferred from the Lagrangian (pl), are collected in Appendix[Aj The vertexes are used in perturbative 
calculations presented in the subsequent sections. 

The Lagrangian (pl) is sometimes [TTJ [T3J E] extended by adding a fundamental N = 2 hypermultiplet to mimic 
a behavior of quarks in QCD plasma. The hypermultiplet is typically massive to study heavy flavors but it can be 
massless as well. We do not consider any extension of the Lagrangian pl) but at the end of Sec. VI we briefly comment 



on a possible structure of self-energies of fields belonging to the fundamental N — 2 hypermultiplet 

III. BASIC PLASMA CHARACTERISTICS 

We start our discussion of SYMP with basic characteristics of the equilibrium plasma. Specifically, we consider the 
energy and partiele densities, Debye mass and plasma parameter of SYMP comparing the quantities to those of QGP. 
For the beginning, however, a few comments are in order. 

In QGP there are several conserved charges: baryon number, electric and color charges, strangeness. The net baryon 
number and electric charge are typically non-zero in QGP produced in relativistic heavy-ion collisions while the total 
strangeness and color charge vanish. Actually, the color charge is usually assumed to vanish not only globally but 
locally as well. It certainly makes sense as the whitening of QGP appears to be the relaxation process of the shortest 
time scale [19]. In SYMP, there are conserved charges carried by fermions and scalars associated with the global 
SU (4) symmetry. One of these charges can be identified with the electric charge to couple N = 4 super Yang-Mills 
to electromagnetic field [3]. In the fortheoming the average SU(4) charges of SYMP are assumed to vanish and so 
are the associated chemical potentials. The constituents of SYMP carry color charges but we further assume that the 
plasma is globally and locally colorless. 

Since there are conserved supercharges in supersymmetric theories, it seems reasonable to consider a statistical 
supersymmctric system with a non-zero expectation value of the supercharge. However, it is not obvious how to 
deal with a partition funetion customary defmed as Tre^^^^^' where /3 = T -1 is the inverse temperature, H is 
the Hamiltonian, Q is the supercharge operator and /i is the associated chemical potential. The problem is caused 
by a fermionic charaeter of the supercharge Q. If \x is simply a number, as, say, the baryon chemical potential, the 
partition funetion even of non-interaeting system does not faetorize into a produet of partition funetions of single 
momentum modes because the supercharges of different modes do not commute with each other. The supercharge is 



TABLE I: Field content of M — 4 super Yang-Mills theory. 

Number of degrees of freedom (N ) 

2 x (Y c 2 - 1) 
6 x (Y c 2 - 1) 
8 x (Y 2 - 1) 



Field's symbol 


Type of the field 


Range of the field's index 


Spin 


.1" 


veetor 


fi, v = 0, 1, 2, 3 


1 


$A 


real (pseudo-)scalar 


A,B = 1,2,3,4,5,6 





Xi 


Majorana spinor 


i, 3 = 1,2,3,4 


1/2 



not an extensive quantity [5D]. There were proposed two ways to resolve the problem. Either the chemical potential 
remains a number but the supercharge is modified by multiplying it by an additional fermionic field c |20| 121 j or the 
chemical potential by itself is a fermionic field [35] . Then, /j,cQ and /j,Q are both bosonic and the partition function 
can be computed in a Standard way. The two formulations, however, are not equivalent to each other. According to 
the former one |20U21j . properties of a supercharged system vary with an expectation value of the supercharge, within 
the latter one |22j . the partition function appears to be effectively independent of Q. Because of the ambiguity, we 
further consider SYMP where the expectation values of all supercharges vanish both globally and locally. 

In view of the above diseussion, SYMP is comparable to QGP where the conserved charges are all zero and so are 
the associated chemical potentials. We adopt the assumption whenever the two plasma systems are compared to each 
other. 

When the chemical potentials are absent, the temperature (T) is the only dimensional parameter, which charaeter- 
izes the equilibrium plasma, and all plasma parameters are expressed through the appropriate powers of T. Taking 
into account that the right numbers of bosonic and fermionic degrees of freedom in SYMP and QGP, the energy 
densities of equilibrium non-interaeting plasmas equal 



_ 30(iV c 2 - 1) 

60 V4(A r < ? - 1) + 7N f N c 



*=^[*,* n \x C .^r.»r T - W 



where the upper expression is for SYMP and the lower one for QGP with Nf light quark flavors. The quark is light 
when its mass is much smaller than the plasma temperature. For N c = N f = 3, the energy density of SYMP is 
approximately 2.5 times bigger than that of QGP at the same temperature. The same holds for the pressure p which, 
obviously, equals e/3. 

The partiele densities in SYMP and QGP are found to be 

2C(3)/ 7(A C 2 -1) \ 3 m 

^ \2(N2-l)+3N f Nj ' {S) 

where £(3) s=s 1.202 is the Ricmann zeta function. For N c = Nf = 3 we have ngYMp/"-QGP ~ 1-3 at the same 
temperature. 

As we show in Sec. |V A| the gluon polarization tensor has exactly the same structure in SYMP and QGP, and 
consequently the Debye mass in SYMP is defined in the same way as in QGP. The masses in both plasmas equal 

2 g 2 f 12N C \ 2 



6 \2N C + N f 

where, as previously, the upper case is for SYMP and the lower one for QGP. For N c = Nf — 3, the ratio of Debye 
masses squared is 2.4 at the same value of gT. The Debye mass determines not only the sereening length r d = 1/wi.d 
but it also gives the plasma frequency u> p = mp/y/S which is the minimal frequency of longitudinal and transverse 
plasma oscillations corresponding to the zero wave veetor. The plasma frequency is also called the gluon thermal 
mass. 

Another important quantity charaeterizing the equilibrium plasma is the so-called plasma parameter A which equals 
the inverse number of partieles in the sphere of radius of the sereening length. When A is decreasing, the behavior of 
plasma is more and more collective while inter-partiele collisions are less and less important. For N c = Nf = 3, we 
have 

As seen, the dynamics of QGP is more collective than that of SYMP. 

The differences of e and n for SYMP and QGP merely reflect the difference in numbers of degrees of freedom in the 
two plasma systems. In case of mp and A it also matters that (anti-)quarks in QGP and fermions in SYMP belong 
to diffcrcnt representations - fundamental and adjoint, respectively - of the SU(iV c ) gauge group. 

IV. DISPERSION EQUATIONS 

Dispersion equations determine dispersion relations of quasi-particle excitations. Below we write down the dispersion 
equations of quasi-gluons, quasi- fermions, and quasi-scalars. 



Gluons 



Sincc the equation of motion of the gluon field A£(k) can be written in the form 



(H) 



where color indices are dropped, H' 11 ' (k) is the retarded polarization tensor and k = (w, k) is the four-momentum, the 
general gluon dispersion equation is 



dct 



fcV" - w - n""(jfe) 



o. 



(12) 



Strictly speaking, one should consider the equation of motion not of the gluon field but of the gluon propagator. 
Then, Eq. (12) determines the poles of the propagator. Due to the transversality of W 1V (k ll H^(k) — 0), which is 
requircd by the gauge covariance, not all components of II'"' are independent from each other, and consequently the 
dispersion equation ( 12 ) can be much simplificd by expressing the polarization tensor through the dielectric tensor 



e %3 (k) which is the 3x3 not 4x4 matrix 



B. Fermions 



The fermion field ipi(k) obeys the equation 

ft -S(fc) t/)(k) = 0, (13) 

where any indices are neglected and S (k) is the retarded fermion self-energy. The dispersion equation thus is 



det 



# -S(fc) 



0. 



Further on we assume that the spinor structure of S(fc) is 

£(fc)=<y*£/.(*). 



(14) 



(15) 



Then, substituting the expression (15) into Eq. (14) and computing the determinant as explained in Appendix 1 of 
[25] . we get 



(fc" -£"(*)) (* M -E M (fc)) 



= 0. 



(16) 



C. Scalars 



The scalar field &A(k) obeys the Klein-Gordon equation 

k 2 + P (k)] $(fc) =0, 



(17) 



where P (k) is the retarded self-energy of scalar field and any indices are dropped. The dispersion equation is 

k 2 + P(k) = 0. (18) 

As seen, the whole dynamical information about plasma medium is contained in the self-energies which are computed 
perturbatively in the next section. 

V. SELF-ENERGIES 



We compute here the self-energies which enter the dispersion equations (12 14 18). The vertexes of N = 4 super 
Yang-Mills, which are used in our perturbative calculations, are listed in Appendix |A[ The plasma is assumed to be 
homogeneous (translationally invariant), locally colorless but the momentum distribution is, in general, different from 
equilibrium one. Therefore, we adopt the Keldysh-Schwinger or real-time formalism which allows one to describe 
both equilibrium and non-equilibrium many-body systems. The free Green's functions, which are labeled with the 
indices +, — , >, <, sym, are collected in Appendix [B| What concerns the Keldysh-Schwinger formalism we follow the 
conventions explained in |23j . The computation is performed within the Hard Loop Approach, see the reviews 
which was generalized to anisotropic systems in [27) . 
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FIG. 1: Contributions to the gluon self-energy. 



Polarization tensor 



The gluon polarization tensor n MI/ can be defined by means of the Dyson-Schwinger equation 

W v (k) = iD^{k) + iD w {k) m pa {k) iV au {k), 



(19) 



where V^ v and D^ v is the interacting and free gluon propagator, respectively. The lowcst order contributions to 
gluon polarization tensor are given by six diagrams shown in Fig. [T] The curly, plain, doted and dashed lines denote, 
respectively, gluon, fermion, ghost, and scalar fields. 

Using the vertexes given in Appendix|A} the contribution to the contour polarization tensor of Keldysh-Schwinger 
formalism, which comes from the fermion loop corresponding to the graph in Fig. [Tk, is immediately written down in 
the coordinate space as 



W n%(x,y) = -ig 2 N c d ab Tr[^S ij (x,y)YS ji (y,x)}. 



(20) 



where the trace is taken over spinor indices. The factor (—1) due to the fermion loop is included and the relation 
facdfbcd = SabN c is used here. 

We are interested in the retarded polarization tensor which is expressed through n > as 



U+(x,y) = e(a; - m )(u>{x,y) -!!<(», y)). 



(21) 



The polarization tensors n > are found from the contour tensor ( 20 ) by locating the argument xq on the upper (lower) 
and j/o on the lower (upper) branch of the contour. Then, one gets 



{ {a) n*(x,y))% = -ig 2 N c 6 ab Tr[ 1 f >Sf j (x,yh"Sf i (y,x)}. 



(22) 



As already mentioned, the system under study is assumed to be translationally invariant. Then, the two-point 
functions as S(x, y) effectively depend on x and y only through x — y. Therefore, we put y = and we write S(x, y) 
as S(x) and S(y,x) as S(-x). Then, Eq. d22h is 



(a) n^(x))^ = --g 2 N c 5 ab Tv[rst 3 {x)^S%{-x)]. 



Since 



S ± (x) = ±Q(±x )(s > (x) - S K (x)) 
the retarded polarization tensor II + (a;) is found as 



(23) 



(24) 



ia) Il + (x))Z = -- 9 2 NA^[rSt 3 {x)YS s ^{-x)+rS s p{x)YSr 3 {-x)], 



(25) 



which in the momentum space reads 



( (a) n + (fc)):: 



,9 



d 4 P 



N c S ab J j^Tr[>fS+(p + k)^S s ^{p) + rS^\p)YSi 3 (p - k)] . 



(26) 



Further on the indcx + is droppcd and II + is d cnotcd as II, as only the retarded pola rization tensor is discussed. 
Substituting the functions S , S sym given by Eqs. (B13 



B15 



B14) into the formula (26 1, one finds 



(a)Kb(k) 



■N c S ab 



d 3 p 2n/(p)-l 



(2tt) 3 



E,, 



(27) 



2p^p v + k^p" + p^k v - g»" (k ■ p) 2p^p u - k»p v - p^k v + g»" (k ■ p) 



{p + k) 2 +isgn((p + k) Q )0^ 



{p -k) 2 -isgn((p- fc) o )0+ 



where p^ = (E pi p) with E p = |p|, the traces of gamma matrices are computed and it is taken into account that 
p 2 = 0. We also note that after performing the integration over p , the momentum p was changed into —p in the 
negative energy contribution. 

In the Hard Loop Approximation, when p 3> k, we have 



1 



1 



2fc 2 



{p + k) 2 + iQ+ (p - k) 2 - i0+ 
1 1 



(k 2 ) 2 - 4(fc • p) 2 - isgn(k ■ p)0+ 
A{k-p) 



1 k 2 

'2{k-p + i0+) 2 ' 

k ■ p 



(p + k) 2 +i0+ (p-k) 2 -i0+ (k 2 ) 2 - 4(fc • p) 2 - isgn(k ■ p)0+ {k-p + i0+) 2 ' 

We note that (p + k) > and (p - k) > for p > k. With the formulas |28||29| ), Eq. ([27]) gives 

d 3 p 2n/(p) - 1 k 2 p^p v - {k»p v +p^k v - g^(k ■ p)) (k ■ p) 



(a) lC(fc) = 4g 2 N c 6 ab 



(2tt) 3 



-E/, 



(fc-p + i0+) 2 



(28) 
(29) 

(30) 



which has the well-known structure of the polarization tensor of gauge bosons in ultrarelativistic QED and QCD 
plasmas. As seen, IL(k) is symmetric with respect to Lorentz indices („jll^^fe) = ( a )Tl vl t(k) and transverse 
^(a)II M r (k) = 0, as rcquired by the gauge invariance. In the vacuum limit, when the fermion distribution func- 
tion n f (p) vanishes, the polarization tensor (30) is still nonzero (actually infmite). As we will see, the vacuum 
contribution to the complete polarization tensor exactly vanishes due to the supersymmetry. 

In analogy to the fermion-loop expression (26), one finds the gluon-loop contribution to the retarded polarization 
tensor shown in Fig. [Tb as 



(&)ICW 



4 



d 4 p f d 4 q 



(2tt) 4 J (2tt) 



D synl (p) \(2ir) 4 6 i4 \k +p - q)M f " y (k,q,p)D + (q) 

+ {27r) 4 S ( - i) {k-p + q)M^(k,-q,-p)D-(q) 



(31) 



where the gluon Green's functions D and D sym are given by Eqs. (BI B4), the combinatorial faetor 1/2 is included 
and 



with 



M^(k, q,p) = F"'(fc, -q,p)T a »(q, -k, -p) 



T^P{k,p,q) = g^{k- P Y + g"P{p-qY + g^{q-kY. 



Within the Hard Loop Approximation the tensor ( 32 ) is computed as 



M" v {k,p ± /c, ±p) « ±2g^{k ■ p) + 10p»p v ± 5(k' J 'p v +p"k u ), 



(32) 
(33) 

(34) 



where we have taken into account that p 2 = 0. 



Substituting the expressions (34) into Eq. (31), using the explicit form of the functions D ± and D sym given by 

(35) 



Eqs. (BI B4), and applying the Hard Loop Approximation (28 29) we get 

,2 



9 



(*«(*:) = ~N c 5 ab 



d 3 p 2n g (p) + 1 5k 2 p^p u - 2g^(k ■ p) 2 - 5(k»p» +p^k l, )(k ■ p) 



(2tt) e 



E, 



i' 



(k -p + iO+y 



The gluon-tadpole contribution to the retarded polarization tensor, which shown in Fig. [Tj:, equals 



(c)KUk) 



9 



d 4 n 

U P T^^Vp n <. 



\ L abccp U UV' 



2 J (2tt) 4 abccp 
where the combinatorial factor 1/2 is included and r^ c p ^ equals 

Kbcd = fa b efec d {9^9 UP ~ <W") + facefedbW V ~ 9^9^) + f ade f eb^ ' 9 P ° ~ ff"V P ). 



With the explicit form of the function D < (p) given by Eq. (B3), the formula (36) provides 

Tl^(b\ 3 n?AT A n^ f ^ 2 %(P) + ! 

(c) U ab W = o9 ^cOabg^ 



(2tt) 3 



E n 



The ghost-loop contribution to the retarded polarization tensor, which is shown in Fig. [TK, equals 



( d )Kb(k) 



-N r S n 



d l p 
(2tt 



4 G sym { P ) {p + ky P »G + {p + k)+p p (p-kyG-{p-k) 



(36) 



(37) 



(38) 



(39) 



where the factor (—1) is included as we deal with the fermion loop. Using the explicit form of the functions G and 
G sym given by Eqs (fB9||B12|), the formula p9l is manipulated to 



{d) Kb(k)=-^N c 5 a b 



9_ 
4 



d 3 p 2n g {p) + 1 feVp" - {W p" +p p k v ){k ■ p) 



(2tt) 3 



E,, 



(k-p + i0+) 2 



(40) 



which holds in the Hard Loop Approximation. 

As already mentioned, the quark-loop contribution to the polarization tensor is symmctric and transverse with 
respect to Lorentz indices. The same holds for the sum of gluon-loop, gluon-tadpole and ghost-loop contributions 
which gives the gluon polarization tensor in pure gluodynamics (QCD with no quarks). The sum of the three 
contributions equals 



(b) + (c)Hd) U ab( k ) = 9 2 N c S ab 



d 3 p 2n g {p) + 1 k 2 p^p u + g^{k ■ p) 2 - (k^p v + p^k v )(k ■ p) 



(2tt) 3 



E,, 



(k-p + io+y 2 



(41) 



To our best knowledge this is the first computation of the QCD polarization tensor in Hard Loop Approximation 
performed in the Keldysh-Schwinger (real time) formalism which explicitly dcmonstrates the transversality of the 
tensor. In Refs. [271 28], where the equilibrium and non-equilibrium anisotropic plasmas were considered, respectively, 
the transversality of LI M1/ (fc) was actually assumed. In case of imaginary time formalism, the computation of the gluon 
polarization tensor in Hard Loop Approximation is the textbook material [291 130) . 

The contribution to the polarization tensor coming from the scalar loop depicted in Fig. [lj is given by 



(e)lC(fc) = -i 3 -8 ab N c S 



AA 



J ^[(2p+kr(2p+kyA + (p+k)A s y m (p) + (2p-kr(2p~kyA s y™(p)A-(p-k)], (42) 



which changes into 



(e) 



n^(fc) = 3g 2 N c 5 ab 



d 3 p 2n s (p) + 1 k 2 pVp v - (p^k v + k^p v )(k ■ p) 



(2tt) 3 



E,, 



{k-p + io+y 2 



(43) 



when the functions A^ and A sym given by Eqs. (B18 B21 ) are used and the Hard Loop Approximation is adopted. 
The contribution to the polarization tensor coming from the scalar tadpole depicted in Fig. [TV is 



(/) 



Kb(k) = -\2ig 2 5 ab N c 5 AA g p » J &A<(p), 



where the combinatorial factor 1/2 is included. With the function A< given by Eq. (B20) we have 

d 3 p 2n s (p) + l 



(/>!!£(*) = 3g 2 N c S ab g^ 



(2tt) s 



E, i 



(44) 



(45) 



a) b) c) 
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FIG. 2: Contributions to the fermion self-energy. 

We get the complete contribution from a scalar field to the polarization tensor by summing up the scalar loop and 
scalar tadpole. Thus, one finds 

TT^m o2ktx f d "P 2n,(p) + 1 feW - {W + WjT - g»" (k ■ p)) (k ■ p) 
(e+f) K b (k) = 3, AU* y ^ — Wp (fc. p + t0+ )2 > ( 46 ) 



which has the structure corresponding to the scalar QED. Thcn, it is not a surprise that the polarization tensor (46) 
is symmetric and transvcrse. 

After summing up all contributions, we get the final expression of gluon polarization tensor 

■»»m_ 2„x / ft/(p)tVl'-(l'p"+lrt'-!"(l-p))('-f) 



HSM - ^A.y ,#X V; + .o4 " (47) 

where 

/(p) = 2n s (p) + Sn/ (p) + 6n s (p) (48) 

is the effective distribution function of plasma constituents. We observe that the coefficients in front of the distribu- 
tions functions %(p), n f (p), n s (p) equal the numbers of degrees of freedom (except colors) of, respectively, gauge 
bosons, fermions and scalars, c/. Table [TJ This is obviously a manifestation of supersymmetry. Another effect of the 
supersymmctry is vanishing of the tensor (H7j) in the vacuum limit when /(p) = 0. Needles to say, the polarization 



tensor (471 is symmetric and transverse in Lorentz indices and thus it is gauge independent. 

In case of QCD plasma, one gets the polarization tensor of the form ( |47[ ) after the vacuum contribution is subtracted. 
For the QGP with the number Nf of massless flavors, the effective distribution function equals 

/qgp(p) = 2n 9 (p) + ^(n g (p) + n g -(p)), (49) 

where n q (p), ng(p) are the distribution functions of quarks and antiquarks which contributc differently to the polar- 
ization tensor than fermions of the JV = 4 super Yang-Mills. This happens because (anti-)quarks of QCD belong to 
the fundamental representation of SU(iV c ) while the fermions belong to the adjoint representation. 

B. Fermion self-energy 

The fermion self-energy E can be defined by means of the Dyson-Schwinger equation 

iS(k) = iS(k) + iS(k) ( - iS(fc)) iS(k), (50) 

where S and S is the interacting and free propagator, respectively. The lowest order contributions to fermion self- 
energy are given by diagrams shown in Fig. [2j The curly, plain, and dashed lines denote, respectively, gluon, fermion, 
and scalar helds. 

The contribution to the fermion self-energy corresponding to the graph depicted in Fig. [2k is given by 



* _2 *r r zin f ^P 

(2tt)< 



(-)S^(*) = \g 2 N c 5 ab 5^ j -[JL [j.S+ip + k)<fir**(p) + llM S s ^(p)^D-(p - k)} . (51) 
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n; 



V: 






d) 



V--- 



FIG. 3: Contributions to the scalar self-energy. 



With the functions D ± , D sym and S ± , S sym given by Eqs. (BI 



B4 



B13 



(o) E£(fc) = 9 2 N c S ab 6^ 



B16), one obtains 
d 3 p n 9 (p) + n/(p) ^ 



(27T) 



-E» 



fc • p + i0+ ' 



(52) 



where the traces over gamma matrices are computed and the Hard Loop Approximation is applied. Eq. ( 52 ) has the 
well-known form of electron self-energy in QED. 

Since there are scalar and pseudoscalar fields X p and Y pi there are two contributions to the fermion self-energy 
corresponding to the graphs depicted in Figs. [2b, [2fc. The first one corresponding to the X p field equals 



d 4 p 



(6) E« (fc) = I y -N c 8 ab a v k ai 3 J ^ [S+ (p + fc)A s >"» + S^{p)A-{p k)] . 



(53) 



Due to the relations (|3|, one hnds that cP ik cP k , = — 3&y. Using the result and substituting the functions S ± , S sym and 



A±, A sym given by Eqs. (B13 



B16 



B18 



B2l|) into Eq. (|53|), one obtains the following result 

i> 



™ m 3 „2 WA hj f d 3 P n/(p) + n s (p) 



(27T) 



-E,, 



k ■ p + i0+ ' 



(54) 



which holds in the Hard Loop Approximation. 

The contribution due to the pseudoscalar field Y p is 



(c) 



E« (*) = z^A^/3^. | ^ [ 75 5+(p + fc) 75 A s ^(p) + 75^ sym (p)7 5 A-b - k)] . 



(55) 



Because /?f fe /3L — ~3$iji 7^75 = — 757/x> an d 75 = 1, we again obtain the result (54 1. 

Summing up all the contributions, we get the final expression for the fermion self-energy 



Ki(fe) = ^N c 5 ab S^ 



d 3 p /(p) j> 

(2tt) 3 £ p fc-p + iOH 



(56) 



which, as the polarization tensor ( 47 ) , depends on the effective distribution function ( 48 ) 



C. Scalar self-energy 

The scalar self-energy P (k) can bc defincd by means of the Dyson-Schwinger equation 

iA(k) = iA{k) + iA(k) iP{k) iA(k), 



(57) 



where A and A is the scalar interacting and free propagator, respectively. The lowest order contributions to the scalar 
self-energy are given by the diagrams shown in Fig. [3] The curly, plain, and dashed lines denote, respectively, gluon, 
fermion, and scalar fields. 

Since there are scalar (X p ) and pseudoscalar (Y p ) fields, we have to consider separately the self-energies of X p and 
Y p . However, one observes that only the coupling of scalars to fermions differs for X p and Y p . The self-interaction and 
the coupling to the gauge field are the same. Therefore, only the fermion-loop contribution to the scalar self-energy, 
which is shown in the diagram Fig. |3k, needs to be computed separately for the X p and Y p fields. 
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In case of the scalar X p field, the diagram Fig. [3k provides 

(a)P^(k) = 4 ^H4 / 7?T4 Tr [^ + (P + k)S^(p) + S s ^(p)S-(p k) 



(58) 



where the symmetry factor 1/2 and the extra minus sign due to the fermi onic character of the loop are included. 
With the explicit form of the functions S , S sym given by Eqs. (B13 B16) and the identity a? -a', = — A5 vq which 
follows from the relations (J3J), onc finds 



(a)P P a !{k) = ~±g 2 N c 5 ab 5™ 



d 3 p 2n/(p) 



(2tt) s 



E„ 



(59) 



The result holds in the Hard Loop Approximation. For the pseudoscalar Y p we obtain the same expression because 
PfiPji =- ~ 4:5 pq , 757^ = —7^75 and j 2 = 1. Therefore, we replace the indices p, q by A, B and we write down the 



result (59 1 as 



(a)Pab(k) = -Ag 2 N c 5 ab 8 



AB 



d 3 p 2n f (p)-l 



(2tt) e 



E n 



(60) 



The contribution represented by the graph depieted in Fig. [Sb equals 



1 



(b)Pab(k) = -1^9 NcSabt 



AB 



d 4 p 



(2tt)' 



O') 



i [(p + 2k) 2 A+(p + k) D sym (p) + (p + k) 2 A s ^(p) D- {p - k)] , (61) 

netions 

Pa b B (k) = ^g 2 N c 6 ab 5' 



which after the substitution of the functions D ± , D sym and A*, A sym in the form pU 

1 



B4 



B18 



B21 ) leads to 



:AB 



d 3 p 4n g (p)-2n s (p) + l 
(2tt) 3 E p 



(62) 



within the Hard Loop Approximation. 

The contributions coming from the gluon tadpole shown in Fig. [3]; and the scalar tadpole from Fig. [3]1 equal, 
respectively, 



(c)Pab (k) = -2g 2 N c 6 ab 6 



AB 



(d)Pab (k) = -5g 2 N c S ab S 



AB 



d 3 p 2n g (p) + 1 

d 3 p 2n s (p) + 1 
(2tt) 3 2E p 



In both cases the symmetry factor 1/2 is included. 

Summing up all contributions we obtain the final formula of scalar self-energy 



P^(fe) = -2g 2 N c 5 ab 6 



AB 



d 3 P /(p) 
(2tt)3 E p ' 



(63) 
(64) 

(65) 



which depends, as II and E, only on the effective distribution funetion (48). 



VI. EFFECTIVE ACTION 



The Hard Loop Approach can be formulated in an elegant and compaet way by introducing the effective aetion 
which was first derived for equilibrium plasmas in [31H33J within the thermal field theory It was also rederived in 
terms of quasiclassical kinetic theory [5H |3S] . Later on a generalization of the aetion to anisotropic systems was given 
in[36j|37]. 

A structure of the effective aetion is constraint by the form of respective self-energies. Since the self-energy of a 
given field is the second functional derivative of the aetion with respect to the field, one writes 



d i yA a Jx)Ii^(x-y)A b l/ (y), 



C 



(*;< 



'(s) 



Af?(z)£S,(x-y)*>) 



4 n) 0*0 = I d*y$ a A (x)P<l»(x-y)<l>%(y) 



(66) 
(67) 
(68) 
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where the self-energies are given by the formulas ( [47] , [56J [65] ) , respectively. The subscript '2' indicates that the above 
effective actions generate only two-point functions. To generate n-point functions these actions need to be modified 
to a gauge invariant form. In the nonAbelian gauge theory studied here, the actions (66 67 68) require a simple 
change - the ordinary derivatives should be replaced by the covariant ones in the final expressions. Rcpcating the 
calculations described in detail in [37J , one finds the Hard Loop effective action of the M = 4 super Yang-Mills as 



C 



HL 



•~*r^+5*i(#*i 






C 



HL 



-HL ' 



where 



r A - 

*-HL — 



*-HL 



s^./tS^W 



f $ 



H L 



g 2 N c 



-2g 2 N c 



(2tt) 3 £ p 

d 3 p /(p) 



^OD^aOD^vOa 



(p- DY I p 



(69) 



(2tt) 3 £1, 
d 3 p /(p 



*?G0 



P'l 



(2tt) 3 £ p 



^(x)$i(x). 



*i 6 (s), 



(70) 
(71) 
(72) 



where /(p) is, as previously, the effective distribution function of plasma constituents. 



The actions (70 71 72) are obtained from the self-energies but the reasoning can be turned around. As argued 



m [32J [33], the actions of gauge bosons (70), fermions (71), and scalars (72) are of unique gauge invariant form. 
Therefore, the structures of hard-loop self-energies of gauge bosons, fermions and scalars are unique. Consequently, 
the self-energies computed in the previous section and those corresponding to the fundamental M = 2 hypermultiplet 
can be inferred from the known QED and QCD results with some help of supersymmetry arguments. However, explicit 
computations, as those presented in Sec.fV] seem to be still needed to determine, at least, numerical coefficients. 



VII. COLLECTIVE MODES 



When the self-energies computed in Sec. [V] are substituted into the dispersion equations presented in Sec. |IV| 
collective modes can be found as solutions of the equations. Below we briefly discuss the gluon, fermion, and scalar 
excitations. 



• The structure of polarization tensor (47) is such as of gluon polarization tensor in QCD plasma. It has also 
analogical form as in both usual and supcrsymmetric QED plasma. Therefore, the spectrum of collective 
excitations of gauge bosons is in all cases the same. In equilibrium plasma we have the longitudinal (plasmon) 
mode and the transverse one which are discussed in e. g. the textbook [53] . When the plasma is out of equilibrium 
there is a whole variety of possible collective excitations. In particular, there are unstable modes, see e. g. the 
review [38] . which exponentially grow in time and strongly influence the system's dynamics. 



• The form of Majorana fermion self-energy (56) happens to be the same as the quark self-energy in QCD plasma. 



It also coincides with the electron self-energy in both non-supersymmetric and supersymmetric QED plasma. 
Therefore, we have identical spectrum of excitations of fermions in all these systems. In equilibrium plasma 
there are two modes of opposite helicity over chirality ratio, see in e. g. the textbook |29) . One mode corresponds 
to the positive energy fermion, another one, sometimes called a plasmino, is a specific medium effect. In non- 
equilibrium plasma the spectrum of fermion collective excitations changes but no unstable modes have been 
found cvcn for an extremely anisotropic momentum distribution |39l I40j . 



The scalar self-energy (65) is independent of momentum, it is negative and real. Therefore, P(k) can be wri tten 
as P(k) 



-m 



cir 



where m c s is the effective scalar mass. Then, the solutions of dispersion equation (17) are 



E,, 



±V~ r 



cir 



+ V Z 



We conclude this section by saying that the gauge boson and fermion excitations of SYMP are the same as in 
ultrarelativistic QED and QCD plasma. The scalar excitations are of the form of free massive relativistic partiele. 
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VIII. COLLISIONAL CHARACTERISTICS 

We consider here characteristics of the Af = 4 super Yang-Mills plasma which are driven by collisions of plasma 
constituents. We start with a review of elementary processes and then we discuss transport coefficients. 

A. Elementary processes 

The elementary processes, which occur at the lowest nontrivial order of the coupling constant g, are binary interac- 
tions, the cross sections of which are proportional to g*. Table|ll|gives the respective matrbc elements squared summed 
over all internal degrees of frccdom of interacting particlcs. G, F, S denote a gluon, fermion and scalar, respectively. 
The matrix elements, which wcre first computed in |10j . are expressed through the Mandelstam invariants s, t and u 
defined in the Standard way. For a process symbolically denoted as 1 + 2 — > 3 + 4, we have 

s = (p 1 +p 2 f, t=(p 1 -p 3 ) 2 , u=(p 1 -p 4 f, (73) 

where Pi,P2,P3,P3 are the four-momenta of particles 1, 2, 3, 4, respectively. For a given process, the differential cross 
section, which is summed over the internal degrees of freedom of final state particles and averaged over the internal 
degrees of freedom of initial state particles, is expressed through the matrix element squared from Table [TT| as 

111 V^ I 71 r|2 (1A\ 

where iVj lof and A^ of are the numbers of internal degrees of freedom of initial state particles given in Table [M The 
collisional processes listed in the Table [TT] determine transport properties of the plasma. 

B. Transport coefficients 

Transport coefficients of weakly coupled QGP, which include baryon and strangeness diffusion, electric charge 
and heat conductivity, shear and bulk viscosity and color conductivity, have been studied in detail, see [HMI] and 
references therein. The shear viscosity of SYMP has been computed in [TD] and the bulk viscosity is identically zero 
because of exact conformality of the system. Other transport coefficients of SYMP have not been studied but one 
expects the coefficients to be qualitatively similar to those of QGP. 

Since the temperature is the only dimensional parameter, which characterizes the equilibrium plasma of massless 
constituents, one finds that, for example, the shear viscosity r\ must be proportional to T 3 and the color conductivity 

TABLE II: Elementary processes in M — 4 super Yang-Mills plasma. 
n Process 

1 GG^GG 

2 GF o GF 

3 GG^ FF 

4 GS-&GS 



5 GG^ S S 

6 GF o SF 

7 GS 1 o FF 

8 FS o FS 

9 SS^ FF 

10 SS^SS 72(^±^ + ^^ + ^^+3) 

11 FF o FF 128f^#^ + 1 ^V- + t - ± ^ L +3) 





i i 

9 4 ^(JV^-l) 


j:\m\ 2 




8(-' 


5 2 +U 2 | U 2 +t 2 
t 2 ' S 2 


, t 2 +s 2 , 
+ n 2 + 


3) 




32 (^+h!_ 


U s \ 
s u J 






32(' 2 +" 2 


u t\ 

t u) 






24(^+h! 


-+1) 






24(^ 


+1) 






-96(7 + 


5 + 1) 






-96(f + 


!+l) 






-96[^ + 3(* + i)+l] 






-96[^+3(i 


- + i) + 1] 
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<7 C to T. It appears that the dominant contributions to both transport coefficients of QGP come from the binary 
collisions driven by a one-gluon exchange which correspond to the matrbc elements squared diverging as t~ 2 for 
t — > 0. The analyses presented in [U] and [U], respectively, show that at the leading order r\ ~ T 3 /g i lng~ 1 and 
a c ~ T/hvg" . The factor 1/lng" 1 appears due to the infrared singularity of the Coulomb-like interaction which is 
regulated by the gluon self-energy. Actually the physics bchind the two formulas is rather different. The viscosity 
is governed by collisions with the momentum transfer of the order of gT while for the color conductivity the softer 
collisions with the momentum transfer of the order g 2 T play a crucial role. 

One expects the same parametric form of rj, a and other transport coefficients in case of SYMP and QGP because, 
similarly to QGP, there are the Coulomb-like binary interactions for every constituent of SYMP, see Table [TTJ. The 
analysis [10] indeed proves that the shear viscosity coefficients of QGP and SYMP differ only by numerical factors 
which mostly reflect different numbers of degrees of freedom in the two plasmas. The viscosity is strongly dominatcd 
by the Coulomb-like interactions, and consequently it does not much matter that the sets of elementary processes in 
the two plasma systems are different. 

In the paper [16] we considered two transport charaeteristies of the TV = 1 QED plasma which are not so constrained 
by dimensional arguments and seemed to strongly depend on elementary process under consideration. Specifically, 
we computed the collisional energy loss and momentum broadening of a partiele traversing the equilibrium plasma. 
The latter quantity determines a magnitude of radiative energy loss of highly energetic partiele in a plasma [45j . The 
dimensional argument does not work here because the two quantities depend not only on the plasma temperature 
but on the energy of test partiele as well. We computed the energy loss and momentum broadening due to the 
processes which, like the Compton scattering on selectrons, are independent of momentum transfer. Such processes are 
qualitatively different from the Coulomb-like interactions dominated by small momentum transfers. We managed to 
obtain the exact formulas of the energy loss and momentum broadening due to the momentum-independent scattering. 
In the limit of high energy of test partiele, which is important in the context of jet suppression phenomenology in 
nucleus-nucleus collisions, the energy loss and momentum broadening appeared to be very similar (at the leading 
order) to those driven by the Coulomb-like interactions. 

The result can be understood as follows. One estimates the energy loss ^ as (AE)/X, where (AE) is the typical 
change of particle's energy in a single collision and A is the particle's mean free path given as A -1 = p a with p ~ T 3 
being the density of scatterers and a denoting the cross seetion. For the differential cross seetion, which is independent 
of momentum transfer, the total cross seetion is a ~ e 4 /s. When a highly energetic partiele with energy E scatters 
on massless plasma partiele, s ~ ET and conscqucntly a ~ e i /(ET). The inverse mean free path is thus estimated 
as A -1 ~ e 4 T 2 /E. When the scattering process is independent of momentum transfer, (AE) is of order E and we 
finally find — 4^ ~ e 4 T 2 . In case of Coulomb interaction we have (AE) ~ — e 2 T, A -1 = e 2 T which provide the same 
estimate of the energy loss. The energy transfer in a single collision is thus much smaller in the Coulomb interaction 
than in the momentum independent scattering but the cross seetion is bigger in the same proportion. Consequently, 
the two interactions corresponding to very different differential cross seetions lead to very similar energy losses. 

We expect an analogous situation in SYMP. There are various elementary process but the energy loss and momentum 
broadening of highly energetic partieles do not much differ from those in QGP. 

IX. CONCLUSIONS 

QCD is obviously rather different than Af = 4 super Yang-Mills theory. Nevertheless QGP and SYMP are surpris- 
ingly similar in the weak coupling regime (at the leading order) . The form of gluon collective excitations is identical 
and the same is true for the fermion (quark) modes. The scalar modes in SYMP are as of massive relativistic partiele. 
The sets of elementary processes are different in QGP and SYMP but the transport coefficients, which are dominated 
by the Coulomb-like interactions, are quite similar. The energy loss and momentum broadening of a highly energetic 
test partiele are also rather similar in the two plasma systems. The differences mostly come from different numbers 
of degrees of freedom in both plasmas which need to be taken into account for a quantitative comparison. 

Acknowledgments 

We are grateful to Simon Caron-Huot for a diseussion on conserved charges in SYMP. This work was partially 
supported by the ESF Human Capital Operational Program and Polish Ministry of Science and Higher Education 
under grants 6/1/8.2.1/POKL/2009 and 667/N-CERN/2010/0, respectively. 



15 

Appendix A: Vertexes of M = 4 Super Yang-Mills 

We collect here the vertex functions which are inferred from the Lagrangian (pl) . Since all fields of the Af = 4 super 
Yang-Mills are, except the ghosts, real, there are no arrows orienting the lines. However, one should remember that 
the momentum of every gluon in the three-gluon coupling is assumed to enter the vertex. In case of the gluon coupling 
to scalars, the momentum of one scalar enters the vertex and the momentum of other one leaves it. 
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Appendix B: Green's functions of N — 4 super Yang-Mills 



We present here the retarded, advanced and unordered free Green's functions of the Keldysh-Schwinger formalism 
which are usually labeled with the indices +,—,>,<, respectively. The system is assumed to be translationally 
invariant and locally colorless. The functions of gluons are given in the Feynman gauge. 



Gluons 



The functions of interest of the free gluon field are of the form 
(D ± ( P )r b = g^6 ab D±(p) = - 2 j" USa , b 

{D>{p)r a l = g fi ' / 5 ab D>{p)=g^S a b^{d(E p -p )[n g (p) + i\ + S(E p + Po ) % (-p)), 

(D^fab = 9^5 ab D<(p)=g^5 ab ^(5(E p -p )n g (p) + 5(E p +p )[n g (-p) + l\), 

(D sym (p)r b = (D > (p)r a b + (D < ( P )r ab ^g^5 ab D s y m (p) 

- g^S ab —[S(E p -p )[2n g (p) + 1] + S(E p + p )[2n g (-p) + l]J, 

where E p = |p| and n g (p) is the distribution function of gluons which are assumed to be unpolarized. The 
is normalizcd in such a way that the gluon density is given as 

p g = 2(N?-l)J^n g (p), 

where the factor of 2 takes into account two gluon spin states. 

One checks that the functions (BI B2 B3) obey the required identity 



D > (p)~D < (p)=D+(p)-D-(p). 



Indeed, the left-hand side of Eq. (B6) equals 



(D > ( P )r b (D<(p)r a , b = ^^ (S(E p - Po ) - S(E p+PQ )) = 2 i ng^6 ab 6(p 2 )(e( Po ) Q(-p )). 



Using the well-known relation 



T - n+ =V-Ti-n5{x), 
x ± iQ + x 



(BI) 
(B2) 
(B3) 

(B4) 
function 

(B5) 

(B6) 

(B7) 

(B8) 



one immediately shows that the right-hand side of Eq. (B6) equals the expression (B7| 



Ghosts 



The functions of the free ghost field are 
GUp) = <UG±(p) 



S a b 



p 2 ± isgn(p o )0+' 
Ga b (p) = 5 ab G > (p) = 5 ab ^(5(E p -p a )[n g (p) + l] + S(E P + p )n g (-p)) , 

G< b (p) = S ab G < (p) = 6 ab ^(5(E p -p Q )n g (p) + 5(E p +p )[n g (-p) + l]y 
G S T(P) = G> b (p) + G<(p)=S ab G*y™(p) 

= 6 ab ^-(5(E p -p )[2n g (p) + 1] + S(E p + Po )[2n g (-p) + l]), 



(B9) 
(B10) 

(BII) 

(B12) 



where n 8 (p) is the distribution function of gluons. The functions (B9 B10 BII I obey the identity analogous to the 
relation (B6|. 
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3. Fermions 



The Green's functions of the free massless fermion field equal 
5±(p) = S^SHp)- M 



p 2 ± isgn(po)0+ ' 
S > (P) = ^— jJ (5(B P - po) [n/(p) - 1] + S(E P +po)n/(-p)J , 



S>(p) = 5. 



^J m (p) = S>( P ) + 5<b) = ^^ m b) 



<%|^ (*(£„ -po)[2n/(p) - 1] + <5(£ p +po)[2n/(-p) - l]). 



(B13) 

(B14) 
(B15) 

(B16) 



where n/ (p) is the distribution function of fermions. The distribution function is normalized in such a way that the 
fermion density equals 



p f = 2(N 2 - 1) 



ri 3 p 

(2tt) 3 



n f (p), 



(B17) 



where the factor of 2 takes into account two spin states of each fermion. The functions (B13 B14 B15) are checked 
to obey the identity S > (p) - S<(p) = S + (p) - S- (p). 



4. Scalars 

The Green's functions of the free massless scalar field are 

A±(p) = Sa^ip) 

A: 



p 2 ±isgn(p )0+' 
^abiP) = 5abA > (p) = -5 ab ^(s(E p - Po )[n s (p) + l] + S(E p + p )n s (-p)^ , 

A< 6 (p) - 6 ab A < (p) = -S ab — [S(E p -p )n 8 {p)+5(E p +p )[n s (-p) + l]j, 
AZ n \p) ^ A> ab (p) + A< b ( P ) = 5 ab A^(p) 

= -6 ab ^-{6(E p -p )[2n s {p) + l] +6{E p +po)[2n.(-p) + l]), 



(B18) 
(B19) 
(B20) 

(B21) 



where n s (p) is the distribution function of scalars. The function is normalized in such a way that the scalar density 
equals 



Ps = {Nj - 1) 



d 3 p 
(2tt) 3 



n s (p) 



(B22) 



The functions (B18 B19 B20) obey the identity such as the relation (B6| 
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